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Abstract 



It is shown, how to generate infinite sequences of differential equa- 
tions of the second order based on some standard equations, using Euler- 
Imshenetsky-Darboux (EID) transformation. For all this, factorizations 
of differential operators and operational identities are used. Some genera- 
lizations of integrable cases of Schrodinger's equations are finded. The ex- 
ample of integrable equation with liouvillian coefficients, that, apparently, 
can not be solved by Singer and Kovacic's algorithm (and its modifica- 
tions), was built. The alogorithm and the program for solving conctructed 
classes of equations, was realized in REDUCE system. Corresponding pro- 
cedure GENERATE is addendum to the ODESOLVE procedure of RE- 
DUCE system. The results of using GENERATE procedure (REDUCE 
3.8) was compared with results of DSolve procedure (Maple 10). Though, 
algorithm, based on EID transformation is not an alternative to the un- 
viersal ones, in the borders of its applicability, its very powerful. 

1 Introduction 

This papaer is dedicated to the classical Euler-Imshenetsky-Darboux transfor- 
mation [1,2,3] and its applications to the solving of differential equation. First 
of all, lets make some general remarks 
Let's consider incomplete equations 



V +a y = 0, 



a {x)eC(I), I={a,b). 



(1.1) 



Any complete equation 



y" + aiy' + a Q y 



0, aiGC 1 ^), a eC(I) 



(1.1' 



can be reduced to (1.1) using the substitution 




Indeed, in this case we shall get incomplete equation 

Y" + A (x)Y = 0, A =ao- \a\- ^a[, 

so, equations of the form (1.1) will be considered hereafter. Euler-Imshenetsky- 
Darboux transformation (EID) is transformation of the form 1 

z = p(x)y' + a(x)y, f3(x), a(x) e C 2 (I). (1.2) 

EID transformation is one of the two most important transformations of the lin- 
ear equations. Another important transformation is Kummer-Liouville's trans- 
formation, that has a form 

y = v(x)z, dt — u(x)dx, 

where v(x) and u(x) — are sufficiently smooth functions (you can read about 
this transformation in detail in [4]). 

Statement of the problem: it requires to bring (1.1) to the form 

z" + b o (x)z = 0, b {x)eC(I), (1.3) 

using reversible transformation (1.2). In other wods, 

• by the giving equations (1.1) and (1.3) find transformation (1.2). 
This problem allows another equivalent formulatings: 

• by the given (1.1) and (1.2) find (1.3); 

• by the given (1.2) and (1.3) find (1.1). 

1.1 Statement of the problem in the matrix form 

Let's pass from the scalar equations (1.1)-(1.3) to the matrix ones. Introduce 
designations 

*-(})• *-(;)■ *-U 0- »-U 

Instead of (1.1) (1.3) we shall get: 

Y' = AY, 

Z = TY, T=( , a R f^V detTVO, 

\ a — pa a + p J 

Z' = BZ. 

It is required to find the transformation(1.6), that reducing (1.5) to (1.7). 
Theorem 1 [5, 6] (see also [4]). Next propositions are equivalent: 

1 Transformation (1.2) is well known as Darboux transformation and also Backlund trans- 
formation. We propose to name it EID, that is tribute to the historicas justice (see [4]). 




(1.5) 
(1.6) 
(1.7) 
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a) Equation (1.1) can be reduced to (1.3) by transformation EID (1.2). 

b) System (1.5) can be reduced to (1.7) by transformation (1.6). 

c) Matrices A, B andT in formulas (1-4), (1-6) are such, that the next condition 

is held 

T' = BT - TA. (1.8) 

d) Functions (3(x) and a(x) satisfies the system of the equations 

a" + {b - a )a - 2a [3' - (3a = 0, (1.9) 
13" + (b - a )p + 2a' = 0. (1.10) 

e) The system of the equations (1.9), (1-10) lets first integral (FI) 

a/3' — (id + a 2 + a /3 2 = K, (1.11) 
or (in equivalent form) 

a/3' + /3a' + a 2 +f3p" + bofi 2 =K, (1.12) 

where K — is arbitrary constant, K ^ 0: otherwise EID transformation 
(see formula (1-6)) will be degenerated. 

f) Equation (1.3) by the transformation 

y = z z (1-13) 

y K K V ' 

can be reduced to (1.1). 

g) System (1.7) by the transformation 



Y = T~ X Z, T- 1 = (1.14) 

\ aofS — a' a_ I 
\ K K ' 

can be reduced to (1.5). 

h) Between differential operators, that characterize direct and reverse EID con- 
version consequently 

P = f3D + a, Q = /3D-a-/3', (1.15) 

there are next commutative relations 

(QP + K)Q = Q(PQ + K), (1.16) 

P(QP + K) = (PQ + K)P, (1.17) 
where K is the first integral of the form (1.11), or (1.12). 
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Theorem 1 connects corresponding results of the works [5,6]. Besides, in 
[6] irreversible transformation (1.2) was also considered. The following proof is 
taken from [4]. 

• a) => b). Immidiate corollary, because of the designations (1.4). 

b) => c). Substituting (1.6) in (1.7), we shall get consecutively T'Y + TY' = 
BTY, T'Y + TAY = BTY, from which come to (1.8). 

c) => d). Let's write equation (1.8) in expanded form 

/?' \ 

a" - f3'a - (3a' a' + (3" ) 

( 1 \ / a (3 \ _ ( a (3 \ ( 1 \ 

\ -a J \ a' - (3a a + (3' J \ a' - (3a a + /?' J \ -b J ' 

from which follows the system of the equation (1.9), (1.10). 

d) => e). Multiply (1.9) by (3, (1.10) by a and subtract first equation from 
second: a/3" — (3a" + 2aa' + 2a,Qf3f3' + a' j3 = 0, from which follows FI of the form 
(1.11). Multiplying (1.10) by (3 and adding with (1.11), we shall get (1.12). 

e) => f). From (1.2) by virtue of (1.1) the next equation follows 

z' = (13' + a)y' + (a' - (3a )y, 

(which coefficients corresponds to the lower row of the matrix T) . 

Solving the system, consisting of (1.2) and finded equation, relative to y, in 
virtue of (1.11) we shall come to (1.13). 

f) => g). Reversing the transformation(1.6), we shall come to (1.14) in virtue 
of the relation 

detT = K; 

g) h). Indeed, (3 2 {D 2 + a a ) = QP + K, (3 2 (D 2 + b ) = PQ + K, where 
K =FI. But in this case commutational correlations (1.16), (1.17) are evident. 

h) a). Condition (1.15) include not only direct, but also reverse transfor- 
mation EID. • 

1.2 On solutions of linear equations 
of second order 

Linear differential equation of the second order (1.1'), that allows factorization 
of the form (D — 012) {D — a\)y = 0, has general solution 

y(x)=eJ aidx ( Cl +c 2 J el {a2 - ai)dx dx). 

In particlular, incomplete equation (1.1) allows the factorization 
(D + a)(D — a)y = and has general solution 

y(x) =J adx (c 1+ C2 J e- 2 I adx dx) 

(more detail about this see in [4]). 
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2 Procedure of «generating» 
of the differential equations 
using EID transfomation 

Later on, instead of (1.2) we shall consider transformation 

z = y'-ay (0 = 1, a -» -a). (2.1) 

From (1.9), (1.10) in virtue of (2.1) follows relations 

a" + (b - a )a + a' = 0, b = a + 2a'. (2.2) 

Equation (1.11) will have a form 

a' + a 2 + a = K (2.3). 

Though, (2.1) is the special case of the transformation (1.2), it is of interest as in 
the examination of the spectral problems both in generation of the remarkable 
sequence of the related equations, generated by the equations (1.1). 

Theorem 2 (Euler-Imshenetsky-Darboux) [4]. Equation (1.1), using trans- 
formation (2.1) generates the following sequence of equations 

Vn + anVn = 0, (2.4) 
n 

a n = a + 2^V s _ 1 , (2.5) 

8=1 

where a s -i satisfies the Riccati's equation 

a 's-l + a l-i + a s-i = A s -i. 

// denote 

n ~ \/ Vs — l 

a s -i = (lny s _i) = , 

Us-l 

where y s -\ is eigenfunction of the equation 

y"-i + - x )Vs-i = 0, 
that corresponding to the eigenvalue A = A s -i, than a n will have a form 



a n = ap + 2 

s=l 



n / ~i \ I 
Vs-\ 



Vs-1 



and solution y n (x) of the equation (2.3) can be represented in one of the following 
three forms: 

1 I-, 1 , - V 

Vn= Il( fl -« ! -i)»= l[(D-^)y = y n - 1 J] [D^ \ y , (2.6) 

„—„ a —„ Us-l V Vs—i J 
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where y = y , a = a, y_i = 1. 

We observe that partially contents of the theorem 3 can be found in [1-3]. 
The method of the proof is taken from [4]. 

• Starting equation (1.1) by the transformation 

jh = (D- a )y =(d-4-) y , 

V yoj 

where ao(x) satisfies the equation 

a' + a 2 + a a = A 
when A = A and j/o is eigenfunction of the equation 

y" + (ao - \)y = 0, (2.7) 
that corresponding to the eigenvalue Ao, is reducing to the form 

y'[ + ( a + 2 Ui = 0. (2.8) 



Next, by the transformation 



where y\ is eigenfunction of the equation 

y»+(a + 2(^j' -\J yi = 0, 
that corresponding to the eigenvalue Ai, (2.7) is reducing to the equation 

y'o \ ' , ofy'i 



Proceeding this, we shall find, that the substitution 



y n =[D-y^)y n ^, (2.9) 



Vn-l, 

where y n -i is eigenfunction of the equation 

Vn-l + (ao + 2 If (4) - A j y n _! = 0, 



fe=0 



corresponding to the eigenvalues A n _i, reduces to the equation 
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that corresponding to (2.3). Substitution (2.9) reduces to the first of the equality 
(2.5). The next representation of the factorization is used in proof of the last 
equality in formula (2.6) (see [4, § 5.5]): 

Y[(D - a s _i) = cxp( / a„_idx) JJ[Z)exp(- / (a s _i - a s - 2 )dx)], a_i = 

s=n J s=n J 



Theorem 2, in fact, gives us an algorithm for generating the sequences of 
the multiplying equations, that can be realized in the interactive mode. The 
algorithm's feature is that it simultaneously construct both equation and form 
of its solution. 



2.1 Example 

We consider the equation 

^' + A 2 j/ = 0. (2.10) 

When A = let y = x. Then y' = 1, y' /y = 1/x, {y' yo)' = -1/x 2 . We shall 
get an equation 

y'; + (\ 2 -^)yi = 0, (2.11) 
it's general solution has a form 

yi(x) = — — ^ (ci cos Ax + c 2 sin Ax). (2-12) 

Accept (2.11) (received from (2.12) when c\ = 0, C2 = 1,A= 1) function 

yi(x) = (D — — ) sinx = cosx — — sinx, 
x x 

as the partial solution. Then 

y[ 1 xsinx f Vi\ 1 x 2 — sin 2 x 



yi x x cosx — sinx' \yi J x 2 (x cos x — sinx) 2 
So, in the result, we shall get an equation 

y'i + (A 2 - , 2( " 2 = S[U ' % ) V 2 = 0. (2.13) 
(xcosx — sini)^ 

General solution of the (2.11) can be represented in one of the next forms: 

i X sin X \ 

2/2(20 = {D H 1 ){D )(ci cos Ax + C2 sin Ax), 

x x cos x — sin x x 



or 



sin x 

J/2 (x) = (D z H : D : ) (c\ cos Ax + c 2 sin Ax) . 

xcosx — smx xcosx — sinx 
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Applying REDUCE 2 , we shall get final expression for 2/2(2;)- 



. . — A 2 x cos Xx cos x + (A 2 — 1) cos Ax sin x — Xx sin Xx sin x 

y 2 {x) = a : h 

xcosx — smx 

—X 2 x cos x sin Xx + (A 2 — 1) sin Ace sin x + Xx cos Xx sin x 

+c 2 : 

x cos x — sin x 

Authors can't solve this equation, using Singer and Kovacic's algorithm [8, 9]. 
DSolve procedure (Maple 10) also can't solve this equation. 

Using theorem 2, one can enlarge the number of examples, but we shall be 
interested in those of the sequences, whose common terms can be counctructed, 
using the identities, that are given below. 



3 Special class 

of the related equations 



Theorem 3 // common term of the equation sequence, generated by the EID 
transformation has a form 



ao + n(n + 1) 



y n = 0, ^ = a , 



then its solution (general or partial) can be represented in the form 



(3.1) 



1 v' 1 

y n = n ( D - k '^)yo = yo(D^) n y , 



or 



y n = n ( D D ^yo = y n o + \y^D) n+1 D-\^ 

k=n ^ VO / 



(3.2) 



(3.3) 



where yo is common or partial solution of the equation (1-1), and yo is partial 
solution of the equation (2.6), corresponding to the eigenvalue Ao. 



• Let's accept, that yu = z7o +1 - Then 



a„ = a + 2 [(Inyo)" + (myi)" + . . . + (In y n -i)"] = a a + n(n + 1) 



And 



y n =[D-n^)...(D-2^) (n ^, . 



yo 



yo 



yo 



2 We have used REDUCE 3.8. 
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This transformation can be rewrited in the another form 

Vn =I[(d- k^) D-^yo = y^\y^D) n +'D- l y . 
Indeed, y n in factorized form can be transformated: 

n< D -*i> -«<"£>"■ 

fc=n 

Let's multiply this identity by D in right 

f[(D- kf) = yZ(D±) n D = y^(L D){ L DT = 

£L„ ya ya yo ya yo 

from where we shall have (3.2), (3.3). • 

Let's find members of the class (3.1), where yo satisfies the equation (2.7) 
where A = Ao. Because factorization D 2 + ao — Ao = (D + ao)(D — ao) takes 
plase, we shall have ao = Ao — a§ — a' . Correspondingly afe = Afe — a\ — a' k . 
And also 

k k 
ak = ao + 2 ^2 a' s _ 1 = A - a% - a' + 2 ^ sa' Q . 

S=l 8=1 

We also consider, that a s _i = sa . After simple calculations 

flfe = A fe - (k + l fal -(k+ l)a' Q . 

As a result we shall have an equation k(k + 2)(a + a 2 ,) = Afe — Ao- When 
Afe = A we can find a' + a^ = 0. This equation has a solution a = l/(x + c), 
c — is an arbitrary constant. Than a (x) = 0, Afe = A = 0, i.e. (2.7) will have 
a form y" — 0. Accepting Afe = (fc + l) 2 , we shall find solutions of the equation 
a 'o + a o = 1) that has a form ao = tanh(.x + c), ao = coth(a; + c). Accepting 
Afe = — (fc + l) 2 , we shall find for equation a' + a 2 , = — 1 solutions of the form 
ao = — tan(x + c),ao = cot (a; + c). In accordance with theorem 3, we shall 
construct some important examples. 

4 Operational identities and some types 
of the Schrodinger's equations 

Let's consider equations, generated by equations 

y'o - Ay = o. (4.1) 

1°. A = 0, y = x, Afe = 0, y k = x k +\ a n = -A - ^±^, i.e. we shall 
have an equation 

tf-fA+n^Wo. (4.2) 
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(*-;) 









By virtue of operational identity 

1 
I 

k—n 

that connects factorization of the differential operators with iteration of the 
same differential operator of the first order, the next identity is also correctly 

{{{D- k -)= X ^\ 1 -DY+\ 

k—n 

One can get it, multiplying previous identity by D in right. Indeed, one can go 
to the required identity, transposing factors in the right part of the identity and 
using the associative rule. So, general (partial) solution of the equation (4.2) 
can be express in the form: 

Vn = x n+1 l-Dj yo, (4.3) 

where yo is general (partial) solution of the equation (4.1). 

General solution of the equation (4.1) and also its partial solutions (with 
concrete values of c\, C2) can be express with exponents, hyperbolic and trigono- 
metric functions: 

yo — ci exp(V\x) + a exp(— VXx), (4.4) 

Vo = c \ cosh(v / Ax) + C2 sinh(v / Ax), A > 0, (4.5) 

yo = c\ cos(V— Ax) + c 2 sin(\/— Xx), A < 0, (4.6) 

where c\ and C2 are arbitrary constants. 

When A = differential operator (DO), corresponding to (4.2), allows fac- 
torizations 



D 2_n(n + 1) = D + n+l\ D _n+1 



And corresponding equation 



_ 2(^+1) = 

X z 

has general solution 

y = c lX n+1 + c 2 x- n (4.8) 

(see § 1.2). However, (4.7) is Euler's equation, and its solution (4.8) can be 
finded, using known Euler's substitution x = e*. 

2°. A = 1, y = coshx, X k = (k + l) 2 , y k = cosh fc+1 x, a n = -A - "^2^ , 
i.e. we shall get an equation 

^-(a-^Wo. (4.9) 
V cosh x J 
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Let's multiply operational identity 



1 / l V 

n(D-k tanh a;) = cosh" x D — — 
V cosh x I 

k=n x ' 



by D in right. Then transpose factors in the right part of the identity, using 
associative rule. We shall get an identity 



o / 1 ^ ™ +1 

TT (D - k tanh x) = cosh" +1 x I —D 

11 V cosh x 

k—n 



(4.9) has a solution 



1 \ " +1 

y n = cosh" +1 x ( —^,D\ y , (4.10) 



where yo can be taken in accordance with formulas (4.4)-(4.6). When A 
(n + l) 2 DO, corresponding to (4.9), allows a factorization 

D 2 - (n + l) 2 + gfti + ij = (D + (n + 1) tanh a:) (D - (n + 1) tanh a;), 
cosh a; 



and differential equation 



y ''-(n + l) 2 y +^±^ y = (4.11) 
cosh x 



has general solution 



,.^ I ( ft+ ,/_* 3J -) ,4,2) 

(see p. 1.2), where integral in formula (4.12) can be expressed in elementary 
functions (see [7], § 1.2): 

dx 



I 



2(n+l) 

2n-2fc+l 



= sech^x + V - 2 " ( r V ( ",o~* + o? ^ cch— 

2n + l| £^(2n-l)(2n-3)...(2n-2fc + l) | 

3°. A = 1, y = sinha;, X k = (k + l) 2 , y fe = sinh fe+1 x, a n = -A - "^VJ , 
i.e. we shall have an equation 

^-fA + =^U = 0. (4-13) 



sinh x 

Let's multiply operational identity 



1 / 1 Y 

IT (D - fccotha;) = sinh"x D—— 
±L V smha;/ 

fc=n x ' 
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by D in right and transpose factors in the right part of the identity, using the 
associative rule. We shall have an identity 

. ^ n+l 

f\ (D - fccothx) = sinh n+1 x — — D 
1 •'■ \ smh x 



(4.13) has a solution 



y n = sinh" +1 x f — —Dj y , (4.14) 

where yo is one of the values (4.4)-(4.6). When A = (n+l) 2 DO, corresponding 
to (4.13), allows factorization 

D 2 - (n + l) 2 - =(£) + („+ l) cothz)(£> - (n + 1) cothz), 

sinh x 

and corresponding differential equation 

y "-(n + l) 2 y-^±H y = (4.15) 
smh a; 



y = sinh n+1 x ( Ci + c 2 / — 
V J sir 



(4.16) 



has general solution 

dx 

sinh 2 (" +1 ) x 

(see § 1.2), where in accordance with [7, p. 110] 

dx 

i 2(n+ ^ x = 

2n+l\ J (2n-l)(2n-3)...(2n-2fc + l) C&Ch 



2n- 



4°. A = -1, y = cosx, A fc = -(k+ l) 2 , = cos fc+1 x, a n = -A 
i.e. we shall have an equation 



\ n(n+l) 
cos 2 x ' 



yZ-(x + ^P±)y n = o. (4.17) 

\ cos^ X J 



Let's multiply operational identity 



l / 1 \ n 

n(D- ktanx) = cos n x D 
V COS X J 

k=n 



by D in right and transpose factors in the right part of identity, using the 
associative rule. We shall have an identity: 



n(D - k tana;) = cos n+1 x ( d] 
\ cos X I 

k=n x ' 



n+l 
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(4.17) has a solution 

y n = cos™ +1 x ( — - — D I y , (4.18) 



COS X 

where yo is one of the values (4.4)-(4.6) When A = — (n + 1) 2 DO, corresponding 
to (4.17), allows factorization 

D 2 + (n + 1) 2 - ^ =(D- (n + l)taxix)(D + (n + 1) tana;), 

cos^ X 

and corresponding differential equation 



y" + (n + l) 2 y-^^y = (4.19) 
cos^ x 



has general solution 



y = co^x[c l+C2 J — (4.20) 



cos 2(n+l) x 

(see § 1.2), where in accordance with [7, p. 149] 

[ dx = 

J cos 2 (" +1 )x 



sinx J 2n+1 >^ 2 k n(n - 1) . . . (n - k + 1) 2n-2fc+i 1 

"2^TT\ SeC * ^ (2n - l)(2n - 3) . . . (2n - 2k + 1) S6C X J ' 

5°. A = -l, yo = sinx, \ k = -(k + l) 2 , y k = S m k+1 x, a n = -X-^^, 
i.e. we shall have an equation: 

y''-(x + ^-)y n = 0. (4.21) 
\ sin x J 

Let's multiply operational identity 



l / 1 \™ 

) T (D - k cot x) = sin™ x D 

±L V sin x I 

k=n 



by Z? in right and transpose factors in the right part of the identity, using 
associative rule. We shall have an identity 



k=n 

Equation (4.21) has general solution 



Y[ {D - k cot x) = sin™+ 1 x ( -^—d) 



n+1 



y n = sin" +1 x ( -^—D ) ,,„. (4.22) 

\S1M 
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where yo is one of the values (4.4)-(4.6). When A = —(n+1) 2 DO, corresponding 
to the (4.21), allows factorization 

D 2 + (n + lf- n( - n +^ = (D + (n+l)cabx)(D-(n + l)cotx), 
sin x 

and corresponding differential equation 



y" + (n + lf V - nv ; ^ tl y = (4.23) 
sin x 



(n+1). 



has general solution 



(see § 1.2), where in accordance with [7, p. 148] 

dx 



sin 2 (" +1 > x 



2n + 1 1 (2n - l)(2n - 3) . . . (2n - 2fc + 1) 

Remark. Equation (2.7) is known as Sturm-Liouville's equation, and also 
linear Schrodinger's equation. Equations (4.2), (4.9), (4.13), (4.17) and (4.21) 
are known as Schrodinger's equations with soluble (integrable) potentials, but 
we have constructed them by the another way. When n = 1 this equations are 
partial cases of the Lame's equation, 

y" - (2p(x) + v)y = D, 

where p(x) is the elliptical function of Weirstrass . This partial cases corre- 
sponds to the elementary values of Weirstrass function. At that parameter A 
for (4.9), (4.13), (4.17) and (4.21) when n = 1 differs from v in Lame's equation. 
In [4, p. 67-70] this was explored, using another methods. 

5 Generalized operational identities 

and corresponding Schrodinger equations 



5.1 Equation 

by virtue of identity 

n( D -*^)-<"+»r M ( = +t J >) (") 
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has general solution 

y n = (ax + b) n+1 ( ——rD ) y , (5.3) 



I \ "+ 1 



ax + b 

where yo is one of the values (4.4)-(4.6). When A — we shall have an equation 

n(n + l)a 2 

Vn (ax + b) 2 y " " U ' {bA) 

that allows factorizations 

d2 _ n(n + l)a 2 = ( D _ na \ ( | na 



(ax + b) 2 \ ax + b J \ ax + b 

(n + l)a 



D , (n+l)a 
ax + b 



D 



ax + b 



and has general solution 

y = c 1 (ax + b) n+1 +c 2 (ax + by n . (5.5) 

5.2. Let 

A = m 2 , y = ae mx + be- mx , \ k = m 2 (k + l) 2 , y k = (ae mx + be~ mx ) k+1 . 
Then 

y' _ m(ae mx - be~ mx ) ( y'o\ 4abm 2 



yo ae mx + be- ,nx ' \y J (ae mx + be~ mx ) 2 ' 

Equation 

„ /, 4abm 2 (n + l)n \ , ,„ . . 

- A - 7 ~, -To )Vn = 0, A ^ (5.6) 

by virtue of identity 
TT ( £> - fcm- p ) = (ae mx + be~ mx )" +1 ( \ D 

(5.7) 



k=n 

has general solution 



= (ae™ + be—) n+1 (-^-^d)^ y , (5.8) 

where yo is one of the values (4.4)-(4.6). 
Let A = m 2 (n + l) 2 . By virtue of 

o 2 o Aabm 2 n(n+1) 

D 2 - m 2 (n + l) 2 + 



(ae™ + be- mx ) 2 



15 



equation 

,2 



n i, , U 2 , Aabm 2 n{n + 1) 
y -m(n + l)y+ {aemx + be _ mx)2 y = Q 



has general solution 



y= (ae mx + be- mx ) n+1 



ci+c 2 J (ae mx + be- mx )' 



2(n+1) dx 



(5.9) 



(5.10) 



Let 

A = m 2 ; j/o = acoshmx + fesinhmx, \k=m 2 (k + l) 2 



jjk = (a cosh ma; + bsmhmx) k+1 



Then 



y a sinh mx + b cosh mi f Vo\' 2 a 2 — &' 



yo m a cosh mx + 6 sinh mx ' \yo/ m (acoshmx + 6 sinh mx) 2 ' 
5.3. Equation 

( n(n + l)m 2 (a 2 - b 2 ) \ 

£-['- -, h —r— h V P" = (5 ' n) 

\ (a cosn mx + b sinn mx) / 

by the virtue of the operational identity 

o 



n/ , asmhmx + frcoshmx \ 
D - km — — : 
\ a cosh mx + o sinh mx I 



( 1 

= (acoshmx + sinh mx)™ I D) (5-12) 

\ a cosh mx + b sinh mx / 

has general solution 

/ 1 \ " +1 

y„ = (acoshmx + 6sinhmx)™ + ( : ; — — D ) yo, (5.13) 

\ a cosh mx + b sinh mx / 

where yo is one of the values (4.4)-(4.6). 

Let A = m 2 (n + l) 2 . By the virtue of the operational identity 

D 2 - m 2 {n+l) 2 ' -(n + l)m 2 (a 2 -b 2 ) 



^ , .a sinh mx + 6 cosh mx 
D + m(n + l) 



(a cosh mx + b sinh mx) 2 

a sinh mx + 6 cosh mx 



a cosh mx + & sinh mx 



D-m(n + 1) 



a cosh mx + & sinh mx 
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Equation 



n(n + l)m 2 (a 2 — b 2 ) 



y" -m 2 (n + l) 2 y + 



(a cosh mi + b sinh mi) 2 



y = 



(5.14) 



has general solution 

y = (a cosh mi + &sinhmi)™ +1 

Let 



Ci + c 2 J (a cosh mx + b sinh mx) 2(%n+1 ^ dx 

(5.15) 



Aq = — m 2 



yo = acosmx + bsinmx, Xk = — m 2 (k + l) 2 , i/fc = (a cos mx + b sin mx) 
Then 



fe+i 



?)q a sin mx + b cos mx { y'o\ 



-m 



a 2 + b 2 



y a cos mx + b sin mx ' \ y 
5.4. Equation 

m 2 (a 2 + 6 2 )n(n + l) 



(a cos mx + b sin mx) 2 ' 



y'n - A + 



(a cos mi + 6 sin mx) 2 
by the virtue of the operational identity 



Vn = 



(5.16) 



n 

k—n 



D — km- 



-a sin mx + b cos mi 



= (a cos mi + b sin mi)™ +1 
has general soltuion 

y n — (a cos mx + b sin mi)" +1 



a cos mi + 6 sin mi 
1 



a cos mi + b sin mi 



-D 



n+l 



a cos mi + b sin mi 



71+1 

D) y , 



(5.17) 



(5.18) 



where yo is one of the values (4.4)— (4.6). 

Let A = — m 2 (n + l) 2 . By the virtue of the factorization 



n(n + l)m 2 (a 2 + b 2 ^ 



D 2 +m 2 (n + l) 2 



(a cos mi + b sin mi) 2 



= D + m(n + l) 
Equation 



-a sin mi + b cos mi 



a cos mi + b sin mi 
y" + m 2 (n + l) 2 y 



£>-m(n + l)- 



-a sin mi + b cos mi 



., . o n(n + l)m 2 (a 2 + fe 2 ) 



a cos mi + b sin mi 
y = (5.19) 



has general solution 

y = (a cos mi + 6sinmi)" +1 



(a cos mi + b sin mi) 2 

i f I | j • N-2(n+l) 

c i + c 2 / (a cos mi + b sin mi) v ' 



fix 



(5.20) 
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6 Algorithm for solving 

of the constructed equations 

We consider following equations: 

equation with rational coefficients of the form: 

(ax + b) z J 

equation with exponential coefficients of the form: 
/ W(n + l)n \ 

Vn \* {aemx + he - mx? )Vn-U, 

equation with hyperbolic coefficients of the form: 

„ / n(n + l)m 2 (a 2 - b 2 ) \ 
y (a cosh mx + o sinh mx) J 

equation with trigonometric coefficients of the form: 

\ (acosmi + Minms)^/ 
Generating equation has a form: 

2/0 - tyo = o. 

Exponential type of the generating equations: 

Vo = c \ exp(v / Ix) + C2 exp(— \flx). 
Hyperbolic type of the generating equation: 

Uo = ci cosh(Vlx) + C2 sinh(v / Ix), I > 0. 
Trigonometric type of the generating equation: 

yo = c\ cos(\f~Ax) + c 2 sin(v / — 1x), I < 0. 

We use the type of the equation (lin, expon, hyp, trig for equations with ra- 
tional, exponential, hyperbolic ant trigonometrical coefficients correspndingly,), 
number n of the equations in the sequence, parameters a, b,m,l, the form of the 
generating equation's solution (expon, hyp and trig for exponential, hyperbolic 
and trigonometric types of solution correspondingly) and arbitrary constants ci 

II C2- 

Input: equation's form and it's solution. 
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Variables 



tp 


The type of the equation 


tpyo 

hline n 


Form of the generating equaion's solution 
The number of the equations in the sequence 


a, b, to, I 


Parameters (look above) 


Cl,C 2 


Arbitrary constants 



1: y p := 

2: if tpyo = expon then 

3: y„ ■= ci exp(y/lx) + c 2 exp(-^x) 

4: else if tpyo = trig then 

5: y n := ci cos(\/^lx) + c 2 sin(\/^7x) 

6: else if tpyo = hyp then 

7: y n :— ci cosh(v / Ix) + c 2 sinh(v^a;) 

8: end if 

9: if o6iii;hh cjiy^aii then 
10: if tp = lin and I 7^ then 

/ 1 

11: Vn := (ax + b) n+1 — D y n 

V ( ax + b ) J 

12: else if tp = lin and I = then 

13: y n := (ax + b) n+1 Cl +c 2 (ax + b)- n 

14: else if tp — expon and I 7^ m 2 (n + l) 2 then 

/ 1 

15: y n := (aexp(mx)-bexp(-mx)) n+1 [ — — -D y n 

\ a exp(TOx) + exp{—mx) ) 

16: else if tp — expon and I = m 2 (n + l) 2 then 

17: y n := (aexp(ma;)+&exp(TOa:))™ +1 (ci+c 2 J(aexp(TOx)+&exp(— mx))~ 2 ( n+1 ' > dx) 

18: else if tp = trig and I ^ —m 2 (n + l) 2 then 

/ x xn+l 

19: y n ■= (acos(mx) + bsm(mx)) n+1 — - — — -D y n 

\ a cos(tox) + 6 sm(TOx) / 

20: else if tp = trig and I — —m 2 (n + l) 2 then 

21: y n = (acos(ma;)+6sin(TOx))™ +1 (ci+c 2 J(acos(mx)+&sin(TOx))~ 2 (™ +1 )<fa:) 

22: else if tp = hyp and I =/= m 2 (n + l) 2 then 

/ 1 y +1 

23: y n = (acosh(TOx)+fesinh(TOx))" +1 — . -D y n 

\ a cosh(TOx) + b sinh(mi) / 

24: else if tp = hyp and I = m 2 (n + l) 2 then 

25: y n = (acosh(mx)+fesinh(TOx))" +1 (ci+c 2 j(acosh(mx)+bsmh(mx))^ 2 ^ n+1 ^dx) 

26: end if 
27: else 

28: if tp = lin then 

29: ty n := (ax + b) n+1 

30: a n :— I + n(n + l)a 2 /(ax + b) 2 

31: a :— a/ (ax + b) 

32: end if 

33: if tp = expon then 

34: ty n := (aexp (mx) + &exp (— mx)) n+1 
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a n := I — Aabn(n + l)m 2 / (aexp (mi) + oexp (-mi)) 2 
a :— Aabm 2 / (a exp (mx) + oexp (-mi)) 2 
if I = m 2 (n + l) 2 then 

y p :— (aexp (mx)+6cxp (— mx) n+1 (ci+c 2 J(aexp(mx)+bex-p(—mx))~ 2 ( n+1 '>dx) 
end if 
end if 

if tp — trig then 

ty n '■= (acos(mx) + 6sin(ma;))™ +1 

a„ := I + n(n + l)m 2 (a 2 + o 2 )/(a cos(t7jx) + osin(mx)) 2 
a := m(— a sin(mx) + ocos(rax))/(acos(m:r) + osin(mx)) 
if / = m 2 (n + l) 2 then 

y p := (acos(mx)+osin(mx))" +1 (ci+c 2 /(a cos(TOx)+6sin(mx))~ 2 (" +1 )dx) 
end if 
end if 

if tp — hyp then 

ty n ■= (acosh(mx) + 6sinh(mx))" +1 

a n := I — n(n + l)m 2 (a 2 — 6 2 )/(acosh(mx) + 6sinh(mx)) 2 
a := m(asmh(mx) + b cosri(mx)) / (a cosh(mx) + osinh(mx)) 
if I = m 2 (n + l) 2 then 

y p := (acosh(mx)+6sinh(mx))" +1 (ci+c 2 J(acosh(mx)+osinh(— mi))" 2 '" +1 )di) 
end if 
end if 

for i := 1 to n do 

tmpi := iayk 
tmp2 := dy n /dx 
y n := tmp 2 - tmp\ 
end for 
end if 



7 Examples 



Now you can see examples of the differential equations and their solutions, 
received with GENERATE procedure (REDUCE 3.8, [10,11]) and DSolve pro- 
cedure (Maple 10, [12]). 

Example 1. Let's consider the following equation: 



y" (*)- (j-2 (cosh(x))- 2 )y(x) = 0. 



Solution, received with DSolve procedure: 

1 



y{x) = 



Vsinh 2x(cosh2x + l) 1 /* 



xci(2cosh2x - 2) 3/4 hyperheom ^ 



i r 

, - cosh 2x + - } -+ 
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V sinh 2x 



c 2 (2cosh2x - 2) 3/4 (2cosh2x + 2) 1/4 x 



x hyperheom ^ 



3 - + 1 -Vi, 3 -- 1 -Vi 

2 2 '2 2 



x(cosh2x + 1) 



1 1\ 
, -cosh2x + - ) x 



As we can see, solution of the comparatively simple equation is expressed in 
the hyperheometric terms. 

And this is solution, received with GENERATE procedure, expressed in the 
explicit form through the elementary functions: 



2/1 = 
Example 



e 2\fix Cl ^ % /y cosn (x) - sinh(x)^ — c 2 ^\/Zcosh(x) + sinh(x)^ 



y" (x)-[l + 2 



e Vlx cos h(x) 



a 2 +b 2 



(acos(x) + b sin (x)) 



y (x) = o. 



DSolve can't finds solution: 



({ 



(- 

\dx 2 



2(a 2 + b 2 ) 



y(x) = DESol ( { ( — Y(x) ) + , v „ , Y{x) , {Y(x)} . 



(a cos a; + b sin x) 2 
GENERATE performs the solution in the elementary functions: 

e 2%//x Cl ^ av / ; C os(a;) — ocos(x) + b\/lsm(x) + asin(x)^ 
e Vlx ( a cog^) _|_ 5 sin(x)) 

c 2 ^— a\/i cos(x) — 6cos(x) — 6\/Isin(x) + asin(x)^ 



2/1 = 

+ 

Example 3: 
Solution, received with DSolve: 



+ 



e^~ lx (acos(x) + fesin(x)) 



/ + 



n(n + l)a z 
(ax + b) 2 



y = 0. 



i \ ^ i ,\-n _^I<22±ii 1 I, i r r. i „V / I(ax + 6)\ 

y(x)~_Cl (ax + b) e « hypergeom M— nj, [— 2 nj, 2 + 



+_C2 (ax + b)" +1 e-^^ 1 hypergeom | [n + 1], [2 + 2 n], 2 ^ ^ + &) 
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Solution, received with GENERATE: 
Example 4- 

v"-(l + ^5-)v = 0. 
V sin x J 

Solution, received with DSolve: 

(-4 cos (2 x) -8 + 2Zcos(2:r) - 2Z)sinh ^ arcsin ^1/2 ^2 cos (2x) + 2 )) 



y(x) = Ci 



-l + cos(2x) 



-3 ^2 cos (2x) + 2 cosh ^a/Z arcsin ^1/2 ^2 cos (2 x) + 2^ ^/-2 cos (2 x) + 2\/Z N 
+ -1 + cos (2 a;) 

/ 3 ((iVl + 1/31 - 2/3) cos(2x) + iVl - 1/31 -4/3) 

+ OS — — x 

I -l + cos(2:r) 

y/-2 cos (2 a; ) + 2 cos ([i\fl- l) arcsin (l/2 ^2 cos (2 x) + 2^ 

-1 + cos (2 a;) + 

-3 ^2 cos (2 x) + 2 sin ^(iv^ - l) arcsin ^1/2 ^2 cos (2 x) + 2^j 

+ -l + cos(2x) X 

((iVl + 1/31- 2/3) cos(2x) - 4/3 -isfl- 1/3 i) \ 
-1 + cos (2x) j 

Solution, received with GENERATE: 

e 2V ~ lx Cl (3cos 2 (a;) - 3\/Icos(V) sin(x) + I s\n 2 (x) + sin 2 (a;)) 

y 2 = n ■ 2, \ + 



+ 



c 2 ^3cos 2 (x) +3Vlcos(x) sin(x) +Z"sin 2 (x) +sin 2 (a;)^ 



8 Conclusion 

Differential equations, that has theoretical or practical significance, can be di- 
vided into equivalence classes by the types of the used transformations. All 
equations, that belongs to the same class are related and stay in relationship 
with standard equation, that has gave them birth. Solutions of all equations, 
belongs to the same class, can be expressed through the solutions of the sample 
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equation. In the terms of diffrential algebra (see [13]) it means, that if con- 
structed equation has liouvillian coefficietns, then it's Picar-Vessiot extension 
belongs to the differential field, generated by the corresponding liouvillian co- 
efficients. Sequential applying of the EID transformation allows to find such 
integrable equations, that can't be solve even in so powerful system as Maple, 
though it uses Singer and Kovacic's algorithm. In our work Schrodinger's equa- 
tions with soluble potentials, that generalaize well-known Schrodinger's equa- 
tions, depending of some parametrs, were constructed, using factorizations of 
differential operators. Solutions (eigenfunctions) of this Schrodinger's equations 
are expressed in terms of equations with constant coefficients. The constructed 
algorithm, that finds solutions of the generated equations, also describes a pro- 
cess of their origin. Though, it realized as GENERATE procedure in REDUCE 
3.8, it can be realized in other systems of the computer algebra. In the bound- 
aries of its applicability, this procedure is not less, than universal procedures. 
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